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The best estimate for the masses of glueballs comes from lattice gauge theory calculations 
which in the quenched approximation show that the lightest glueball has J^*^ = 0+"*" and 
that its mass should be in the range 1.45 — 1.75 GeV. While the lattice remains immature for 
predicting glueball decays, Amsler and Close have noted that in lattice inspired models, such 
as the flux tube [0], glueballs will mix strongly with nearby qq states with the same 
This will lead to three isoscalar states of the same J^*-^ with predictable patterns of decay 
branching ratios 0, ^, ^. Such mixing ideas have been applied recently to the three states in 
the glueball mass region - the /o(1370), /o(1500) and /o(1710) 0, |, 0, |, |]. 

While these papers at first sight differ in detail, nonetheless their conclusions share some 
common robust features. The flavour content of the states have the nn and ss in phase (SU(3) 
singlet tendency) for the /o(1370) and /o(1710), and out of phase (octet tendency) for the 
/o(1500) (that such a pattern is natural is discussed in ref. [|13])- In general these mixings 
will negate the naive folklore that glueball decay branching ratios would be independent of the 
quark composition of the final state mesons after taking into account phase space effects, so 
called "flavour blind decays". 

Recently ref. ^ has used some of the observed branching ratios as input to constrain the mix- 
ing pattern. The results here too agree with the generic structure found in refs. 0, |, ^ |10|. 
Most recently, and subsequent to the above work, the WA102 collaboration has now pub- 
lished for the first time in a single experiment, a complete data set for the decay branching 
ratios of the /o(1370), /o(1500) and /o(1710) to all pseudoscalar meson pairs. These data will 
be our point of departure. Using methods similar to those proposed in ref. ^ we investigate 
the implications of the WA102 data for the glueball-quarkonia content of the /o(1370), /o(1500) 
and /o(1710). This moves the debate forwards in the following ways: 



• It highlights the sensitivity of the mixings to the input data. 

• It exposes some assumptions, both explicit and implicit, in the analysis of ref. that 
can be improved upon. 

• It allows for the direct decay of glueballs into t] and r^', which were not manifested in lead- 
ing order in ref. |§] even though there are reasons to suspect that they could be important 

B0- 



In order to unfold the production kinematics we use the invariant decay couplings (7jj) for 
the observed decays, namely we express the partial width (Fj^) as 

T.,=^l\F,m'Sp{q) (1) 

where Sp{q) denotes the phase space and Fij{q) are form factors appropriate to exclusive two 
body decays. Here we have followed ref. and have chosen 

\F,M\' = Q''exp{-q'/8p') (2) 

where I is the angular momentum of the final state with daughter momenta q and we have used 
f5 = 0.5 GeV/c 0. The /o(1500) lies very near to threshold in the rjr]' decay mode, therefore 
we have used an average value of q (190 MeV/c) derived from a fit to the rjr]' mass spectrum. 
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The branching ratios measured by the WA102 experiment for the /o(1370), /o(1500) and 
/o(1710) are given in table |I]. The invariant decay couphngs {'jij) are related to the relevant 
decay amplitudes by 

7j = Q,|M,,f (3) 

where Cij is a weighting factor arising from the sum over the various charge combinations, 
namely 4 for KK, 3 for tttt, 2 for rjrj' and 1 for 7777 for isoscalar decays. If in the decay of 
some state the ratios of the decay amplitudes squared (|Mjjp) for the rirj/nn and rjri/KK are 
simultaneously greater than unity, then this state cannot be a quarkonium decay (see fig. 3 of 
ref. [0] and also ref. 0). Table |^ gives these ratios for the /o(1370), /o(1500) and /o(1710), as 
abstracted using eqs.(|I|), (0) and (0). As can be seen the ratios for the /o(1710) argue either 
for non-gg content in this meson or for some further dynamical suppression of the tttt mode, 
say, as may occur for special values of parameters in some specific models . 



In the \G) = \gg), \S) = \ss), \N) = \uu + dd)/^/2 basis, the mass matrix describing the 
mixing of a glueball and quarkonia can be written as follows M: 



M 




(4) 



where / = {G\M\S) = {G\M\N) / \/2 represents the flavour independent mixing strength be- 
tween the glueball and quarkonia states. Mg, Ms and represent the masses of the bare 
states \G), \S) and |A^), respectively. Following refs. 0, ^ |^ we assume the mixing is strongest 
between the glueball and nearest qq neighbours. With the lattice (in the quenched approxima- 
tion) predicting the glueball mass to be in the 1.45 — 1.75 GeV region, this naturally led these 
papers to focus on the physical states |/o(1710)), |/o(1500)) and |/o(1370)) as the eigenstates of 
M with the eigenvalues of Mi, M2 and M3, respectively. (An alternative picture could involve 
the states /o(1500), /o(1710), /o(2000); we do not discuss this in the present paper). The three 
physical states can be read as 
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(6) 



Ref.|0| considered the following three hadronic decay paths for the /o(1370), /o(1500) and 
/o(1710): 

(i) the direct coupling of the quarkonia component of the three states to the final pseu- 
doscalar mesons (-P-P) (fig. |l]a), 

(ii) the coupling through two intermediate gluons, nn{ss) ^ gg ^ ss{nn) PP (fig. [l|b), 
with ri representing the ratio of the effective coupling strength of this mode to that of the 
mode (i); 
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(iii) the flavour independent coupling of the glueball component gg qiqi with subsequent 
decay q-iql — PP (fig. |l|c); with r2 representing the ratio of this mode to (i) 

We propose that this is inconsistent. Specifically, the modes (ii) and (iii) as described above 
are what have already been accounted for in generating the mixed states and so ri and r2 should 
be set to zero. This may be seen by comparing the definitions of eq.(^ with the corresponding 
expressions in perturbation theory (as e.g. eqs. (27-31) in ref extended to second order). For 
example, eq.(27) of ref. (where Cq denotes the normalisation factor) 



(r\-U(r\ , \ss){s-s\M\G) , \nn){nn\M\G) 

{Co) fo{G)-\G)+ + (7) 



may be written in the form of eq.(6), with Mq = Mi, and N = nn, S = ss 



{MG-M^^){MG-M,,)fo{G) = {Mi-M^){Mi-Ms)Ci\G) + {Mi-MN)fCi\ss) 

+{Mi - Ms)V2fGi\ nn) (8) 

This shows how the coefficients yi, z\ are equivalent to the perturbation which is in turn driven 
by fig. |I]c). Thus it is double counting to invoke this same figure with strength "r2" to describe 
gg decays via qq intermediate states, having already used it to compute the mixing of those 
same qq in the Fock state. Similar remarks apply to the ss — > nn mixing in second order 
perturbation theory and fig. [l|b defined as r^. 

However, there are additional pathways that have not been allowed for in ref. p|. First 
there is the role of gg qqqq as in fig. These may have important resonant contributions 
in the region below 1 GeV but are expected to be primarily a continuum in the 1.5 GeV region 
of interest here |T^; we shall approximate them by assuming flavour independent couplings. 



(In a more sophisticated analysis one could incorporate threshold effects in the PP channels 
that overlap the qqqq configurations |T3; we do not discuss this further in this first look). The 



resulting amplitudes can be obtained from eqs.(A4) of ref. and have the same structure as 
those of (iii) above. Hence a non-zero r2 is restored, though its interpretation differs from ref 



Finally, following ref 0, 12], we allow for 

(iv) the direct coupling of the glue in the initial state to isoscalar mesons (i.e. rjrj and rjri' 
decays). As in ref [||, we assume chiral symmetry such that the coupling to the s's content of 
the 7], 7]' dominates and allow to be the ratio of mode (iv) to (i). 

The three decay diagrams considered are shown in fig. |^-c. Performing an elementary 
SU(3) calculation gives the reduced partial widths in table ^ where a = (cos 0- V2sin 0)/\/6, 
(3 = (sin + a/2 cos (p) / ^/6, and (p is the mixing angle of r] and r]f. The relevant expressions 
follow from appendix A in ref. with p = R = 1 in the case of fiavour independence of the 
direct couplings. The predicted branching ratios can then be calculated using eqs.(|l]) and (0). 

We then perform a fit based on the branching ratios given in table [l|, where we have 
required that the matrix U in eq. (13) is unitary, which applies an additional 6 constraints to the 
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fit. As input we use the masses of tlie /o(1500) and /o(1710). In tliis way eiglit parameters, 
Mtv, Ms, M3, /, r2, and are determined from tlie fit. Tlie parameters determined 
from tlie solution with the lowest ^ire presented in table ^ and the fitted branching ratios 
together with the contributions of each are given in table 0. 

The physical states |/o(1710)), |/o(1500)) and |/o(1370)) are found to be 

|/o(1710)) = 0.39|G) + 0.91|5) + 0.14|Ar), (9) 

|/o(1500)) = -0.691^) + 0.37|^) - 0.62|A^), (10) 
|/o(1370)) = 0.60|G) - 0.13|5) - 0.79|A^). (11) 

It is interesting and non-trivial that the pattern of decays determines flavour mixing angles 
such that a state having an "octet" tendency is sandwiched between two states that have a 
"singlet" tendency. As noted above and elsewhere 0] this is a potential signal for G mixing 



with a qq nonet. The output masses for Mtv and Ms are consistent with the K*{1A?>Q) being 
in the nonet and with the glueball mass being at the lower end of the quenched lattice range 
(see also |TT|). The mixing strength also is in accord with lattice estimates [|^. 



The elements in eqs. (P-pT]) form the matrix U as defined in eqs.(|,|]). We have calculated 
the error on each element taking into account the correlated errors on their constituents which 
gives 

/ 0.14 

Af/ = 




0.07 0.06 0.08 I . (12) 
V 0.08 

This shows that the "singlet-octet-singlet" phase pattern is robust. The most sensitive probe of 
flavours and phases is in 77 couplings. In the spirit of ref. ignoring mass-dependent effects, 
the above imply 



r(/i(1710) 77) : r(/i(1500) ^ 77) : r(/i(1370) ^ 77) = 
(5zi + V2yif : (5^2 + V2y2f : (5^3 + ^Vz? = 3.8 : 6.8 : 16.6. (13) 

The 77 width of /o(1500) exceeding that of /o(1710) arises because the glueball is nearer 
to the N than the S. Contrast previous works where the G was nearer to (or even above) the 
S, in which case the /o(1500) has the smallest 77 coupling of the three states 0. This shows 
how these 77 couplings have the potential to pin down the input pattern. 

An interesting feature is the small value of the pseudoscalar mixing angle (0); it is interesting 



that this value agrees with recent work that has (f){ri) 7^ 4>{ti') [18 1). We have checked that our 
results are not sensitive to allowing these angles to be independent. If instead we fix the value 
of (j) to -19° degrees, the of the fit increases from 3.0 to 7.7 and the results are given in 
tables ID and ^ respectively. As can be seen the parameters of the fit are not very much affected. 

We have also tried setting to zero: the of the fit increases to 13.9 (see tables |I] and ^ 
and once again the parameters are not much altered. 
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Other authors have claimed that Mq > Ms > M^r 19, 20|. This scenario is disfavoured 



as, if in the fit we require Mq > Ms > M^r, the increases to 57. In any event, we are cautious 
about such claims |T9|, ^ as they are likely to be significantly distorted by the presence of a 



higher, nearby, excited nn state (N*) such that Mat* > Mg > Ms: the philosophy of dominant 
mixing with the nearest neighbours would then lead again to the "singlet - octet - singlet" 
scenario that we have found above. We defer detailed discussion to a more complete report. 

Our preferred solution has two further implications for the production of these states in 
pp annihilations, in central pp collisions and in radiative J/ip decays that are in accord with 
data. These are interesting in that they are consequences of the output and were not used as 
constraints. 

The production of the /o states in pp ^ n + fo is expected to be dominantly through 
the nn components of the /o state, possibly through gg, but not prominently through the ss 
components. (The possible presence of hidden ss at threshold, noted by |^ is in general 
swamped by the above, and in any event appears unimportant in flight). The above mixing 
pattern implies that 



aipp ^71 + /o(1710)) < a{pp -> vr + /o(1370)) ~ a{pp ^ n + /o(1500)) (14) 
Experimentally p2i the relative production rates are. 



pp-^n + /o(1370) : tt + /o(1500)) -1:1. (15) 

and there is no evidence for the production of the /o(1710). This would be natural if the 
production were via the nn component. 

For central production, the cross sections of well established quarkonia in WA102 suggest 



that the production of ss is strongly suppressed relative to nn. The relative cross sections 
for the three states of interest here are 



PP^PP+ (/o(1710) : /o(1500) : /o(1370)) ~ 0.14 : 1.7 : 1. (16) 
This would be natural if the production were via the nn and gg components. 

In addition, the WA102 collaboration has studied the production of these states as a function 
of the azimuthal angle 0, which is defined as the angle between the pt vectors of the two 
outgoing protons. An important qualitative characteristic of these data is that the /o(1710) 
and /o(1500) peak as — >• whereas the /o(1370) is more peaked as — > 180 [^]. If the 
gg and nn components are produced coherently a.s (p but out of phase as (p —>■ 180, then 
this pattern of (p dependence and relative production rates would follow; however, the relative 
coherence of gg and nn begs a dynamical explanation. 

In J/ip radiative decays, the absolute rates depend sensitively on the phases and relative 
strengths of the G relative to the qq component, as well as the relative phase of nn and ss within 
the latter. The general pattern though is clear. Following the discussion in ref. we expect 
that the coupling to G will be large; coupling to qq with "octet tendency" will be suppressed; 
coupling to qq with "singlet tendency" will be intermediate. Hence the rate for /o(1370) will be 
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smallest as the G interferes destructively against the qq with "singlet tendency" . Conversely, 
the /o(1710) is enhanced by their constructive interference. The /o(1500) contains qq with 
"octet tendency" and its production will be driven dominantly by its G content. If the G mass 
is nearer to the than to the S, as our results suggest, the G component in /o(1500) is large 
and cause the J/ip — > 7/o(1500) rate to be comparable to J/iIj — > 7/o(1710). 

In ref. [||], the branching ratio of BR{J /i) 7/0) (/o ^ tttt + KK) for the /o(1500) and 
/o(1710) is presented. Using the WA102 measured branching fractions |TI| for these resonances 
and assuming that all major decay modes have been observed, the total relative production 
rates in radiative J /ip decays can be calculated to be: 

W ^ /o(1500) : J/i) /o(1710) = 1.0 : 1.1 ± 0.4 (17) 

which is consistent with the prediction above based on our mixed state solution. 

In these mixed state solutions, both the /o(1500) and /o(1710) have nn and ss contributions 
and so it would be expected that both would be produced in n^p and K~p interactions. The 
/o(1500) has clearly been observed in 7r~p interactions: it was first observed in the 7777 final 
state, although at that time it was referred to as the (^(1590) p6|. There is also evidence 
for the production of the /o(1500) in K~p K^K^A [^, The signal is much weaker 
compared to the well known ss state the /2(1525), as expected with our mixings in eq.(p!OD and 
the suppressed KK decay associated with the destructive nn — ss phase in the wavefunction. 

There is evidence for the /o(1710) in the reaction 7r~p — > KgKgn, originally called the 
5*'(1720) [H], One of the longstanding problems of the /o(1710) is that inspite of its 

dominant KK decay mode it was not observed in K~p experiments ES, However, these 



concerns were based on the fact that initially the /o(1710) had J = 2. In fact, in ref. ||3^ it 
was demonstrated that if the /o(1710) had J = 0, as it has now been found to have, then the 
contribution in 7r~p and K~p are compatible. One word of caution should be given here: the 



analysis in ref. ||32| was performed with a /o(1400) not a /o(1500) as we today know to be the 
case. As a further test of our solution, it would be nice to see the analysis of ref. repeated 
with the mass and width of the /o(1500) and the decay parameters of the /o(1710) determined 
by the WA102 experiment. 



In summary, based on the hypothesis that the scalar glueball mixes with the nearby qq 
nonet states, we have determined the flavour content of the /o(1370), /o(1500) and /o(1710) by 
studying their decays into all pseudoscalar meson pairs. The solution we have found is also 
compatible with the relative production strengths of the /o(1370), /o(1500) and /o(1710) in pp 
central production, pp annihilations and J/ip radiative decays. 
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Table 1: The measured and predicted branching ratios with the individual contributions 
coming from the fits. 





Measured 


All free 


(p - 


1 o° 
-iy 


^3 = 


n 
U 




Branching 


Fitted 




Fitted 




Fitted 


x' 




ratio 














/o(1370)^7r7r 
fo(1370)^KK 


2.17 ± 0.9 


2.13 


0.001 


2.1 


0.004 


2.25 


0.007 


/o(1370)^W 
/o(1370)->KJ£" 


0.35 ± 0.21 


0.41 


0.1 


0.01 


2.6 


0.01 


2.6 


/o(1500)^7r7r 
/o(1500)^W 


5.5 ± 0.84 


5.60 


0.01 


5.6 


0.02 


6.20 


0.69 


/o(1500)^KX 
/o(1500)^7r7r 


0.32 ± 0.07 


n 37 




0.32 


0.001 


0.35 


0.22 


fni 1500) ^rm' 
/o(1500)^j?»7 


0.52 ± 0.16 


0.60 


0.23 


0.5 


0.01 


0.20 


3.9 


/o(1710)-»7r7r 
fo{1710)^KK 


0.20 ± 0.03 


0.19 


0.05 


0.20 


0.002 


0.19 


0.08 


/o(1710)^W 
M1710)~^KK 


0.48 ± 0.14 


0.29 


1.9 


0.17 


4.9 


0.13 


6.1 


/o(1710)^W 
/o(1710)-^W 


< 0.05 (90%c/) 


0.034 


0.27 


0.04 


0.05 


0.06 


0.05 
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Table 2: The ratio of decay amplitudes squared. 





r]r]/7T7i 


rjrj/KK 


/o(1370) 
/o(1500) 
/o(1710) 


0.74 ± 0.51 
0.68 ± 0.11 
7.9 ± 2.4 


1.64 ± 0.96 
2.42 ± 0.64 
1.96 ± 0.64 



Table 3: The theoretical reduced partial widths. 





vi) 


2[2aP{zi-V2yi) + 2aPxirs]'' 






[2a^ Zi + 2V2p\ + r2Xi + 2l3^Xir^f 




— >■ tttt) 




i\fi - 


KK) 


A{\{zi + ^/2yi)+r2X^]'' 
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Table 4: The solutions for the minimum . 



Parameters 


All Free 


(f) = -19° 


r3 = 


9 


3.0 


7.7 


13.9 


Mg (MeV) 


1440 ± 16 


1433 ± 19 


1437 ± 15 


Ms (McV) 


1672 ± 9 


1668 ± 8 


1672 ± 13 


Afjv (MeV) 


1354 ± 28 


1366 ± 25 


1368 ± 29 


M3 (MeV) 


1256 ± 31 


1251 ± 18 


1264 ± 14 


/ (MeV) 


91 ± 11 


95 ± 13 


90 ± 11 


<t> (Deg) 


-5 ± 4 


-19 


-25 ± 4 


r2 


1.02 d= 0.14 


0.92 ± 0.18 


0.95 ± 0.12 


rs 


1.04 ± 0.24 


0.77± 0.26 
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Figures 



Figure 1: Possible Decays to Pseudoscalar meson pairs (PP). a) The direct coupling of the 
qq to the PP pair, b) the coupling of the qq to PP via intermediate gluons, c) the coupling of 
the glueball component to PP and d) the direct coupling of gluons to isoscalar mesons. 



Figure 2: The Decays to Pseudoscalar meson pairs (PP) considered in this analysis, a) The 
couphng of the qq to the PP pair, b) the couphng of the glueball component to PP and c) the 
direct coupling of gluons to the gluonic component of the final state mesons. 
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Figure 2 
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